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Abstract 

This paper constructs tridiagonal random matrix models for general {(5 > 0) /3-Hermite 
(Gaussian) and /3-Laguerre (Wishart) ensembles. These generalize the well-known Gaussian 
and Wishart models for (3 = 1,2,4. Furthermore, in the cases of the /3-Laguerre ensembles, 
we eliminate the exponent quantization present in the previously known models. 

We further discuss applications for the new matrix models, and present some open prob- 
lems. 
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1 Introduction 



1.1 Overview 



Classical Random Matrix Theory focuses on the random matrix models in the following 3x3 
table: 





Real, (3 = 1 


Complex, (3 = 2 


Quaternion, (3 = 4 


Hermite 


GOE 


GUE 


GSE 


Laguerre 


Real Wishart 


Complex Wishart 


(Quaternion Wishart) 


Jacobi 


Real MANOVA 


Complex MANOVA 


(Quaternion MANOVA) 



The two entries in parentheses (in the third column) correspond to less-studied random 
matrix models; the others are mainstream and have been extensively researched and publi- 
cized. The three columns correspond to Dyson's "threefold way" (3 = 1,2, and 4; the three 
rows correspond to the weight function associated to the random matrix model. Other weight 
functions have also been considered (for example, the uniform weight on the unit circle corre- 
sponds to the circular ensembles). 



Zirnbauer 33] and Ivanov |12j produced a more general taxonomy of random matrix mod- 
els. Their characterizations ("tenfold", respectively "twelvefold") are based on symmetric 
spaces, and include Hermite, Laguerre, and Jacobi cases, and also the circular ensembles 
(each of their models can be associated with (3 = 1,2 or 4). 

We propose a random matrix program of study that would generalize (3 beyond the above 
threefold way, thus generalizing the 3x3 cartesian product to 3 x oo, making the leap from 
discrete characterizations to continuous ones. A step in this direction has been initiated by 
Forrester [Q], [10], who studied the /^-ensembles in connection with multivariate orthogonal 
polynomials and Calogero-Sutherland-type quantum systems. Furthermore, in the case of 
the classical Laguerre and Jacobi models, our program goes beyond the quantized exponents 
forced by the classical models, and proposes continuous ones. 

For the benefit of the reader we have expanded the 3x3 table with detailed information 
in Figure [l]. 



1.2 Background 

The Gaussian (or Hermite) ensembles arise in physics, and are identified by Dyson Q by the 
group over which they are invariant: Gaussian Orthogonal or for short GOE (with real entries), 
Gaussian Unitary or GUE (with complex entries), and Gaussian Symplectic or GSE (with 
quaternion entries). The Wishart ensembles arise in statistics, and the three corresponding 
models could be named Wishart real, Wishart complex, and Wishart quaternion. 
The three Gaussian Ensembles have joint eigenvalue probability density function 

HERMITE: f p (\) = 4]J |A; - A/e"^ A ?/ 2 , (1) 

i<j 

with 3 = 1 corresponding to the reals, j3 = 2 to the complexes, /3 = 4 to the quaternions, and 



2 



with 



The best references are Mehta [18] and the original paper by Dyson 

Similarly, the Wishart (or Laguerre) models have joint eigenvalue p.d.f. 

LAGUERRE: f p {\) = c P £ a J] \\ - A,f J] \ a ~ v e~ ^ , (3) 



i<j 



with a = and p = 1 + ^(m — 1). Again, (3 = 1 for the reals, (3 = 2 for the complexes, and 
(3 = A for the quaternions. The constant 



r(i + f) 

1 " ll r (l + fi)r(a-f(m-i))) 



/3,a 



TT ililLlj (4) 



Good references are [£l]], j8|, and [|L|, and for (3 = 4, 17]. 

To complete the triad of classical orthogonal polynomials, we will mention the /3-MANOVA 
ensembles, which are associated to the multivariate analysis of variance (MANOVA) model. 
They are better known in the literature as the Jacobi ensembles, with joint eigenvalue p.d.f. 



JACOBI: f^X) = c^ 2 ]J \X t - A/ J] (1 - KT 2 ~ P , (5) 

i<j j=l 

with oi = %nx, a>2 = §^2, and p = 1 + 4(m — 1). As usual, (3 = 1 for real and (3 = 2 for 
complex; also 



ftax,a 2 = TT £(1 + f )r( Ql + a 2 - j (m - j)) 

J tAni + lmai-§(rn-j))T(a 2 -§(m-j))' 



c 



The MANOVA real and complex cases {(3 = 1 and 2) have been studied by statisticians (see 

0). ^ ^ _ ^ ^ ' 

Though "Gaussian" , "Wishart" , and "MANOVA" are the traditional names for the three 
types of /^-ensembles, we prefer the sometimes used and technically more informative names 
"Hermite", "Laguerre", and "Jacobi" ensembles. These technical names reflect the fact that 
the p.d.f.'s for the ensembles correspond to the p.d.f.'s etr(— A 2 /2), det(^4) a_p etr(— A/2), and 
det(j4) ai-p det(/ — A) a2 ~ p over their respective spaces of matrices. In turn, these functions 
correspond to three sets of orthogonal polynomials (Hermite, Laguerre, Jacobi). Throughout 
this paper, we will use the term "general /3-Hermite, -Laguerre, -Jacobi ensembles" for general 
(3 in the p.d.f.s ©, (|), ©. 

Though it was believed that no other choice of (3 would correspond to a matrix model 
constructed with entries from a classical distribution, there have been studies of general (3- 
Hermite ensembles as theoretical eigenvalue distributions. They turn out to have important 
applications in lattice gas theory (see |1(J, p). 

The general (3 ensembles appear to be connected to a broad spectrum of mathematics 
and physics, among which we list lattice gas theory, quantum mechanics, and Selberg-type 
integrals. Also, the (3 ensembles are connected to the theory of Jack polynomials (with the 
correspondence a = \ where a is the Jack parameter), which are currently objects of intensive 
research (see pi], pj, p§). 
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1.3 Our Results 



Dyson's original threefold way is a byproduct of the invariance assumptions as in the "In- 
variance" column of Figure |]. By necessity, any invariant distribution is generically dense. 
Further, the invariance approach forces the consideration of the complex and quaternion di- 
vision algebras. 

In this paper, we drop the invariance requirement. What we gain are "sparse" models (with 
only O(n) non-zero parameters) over the reals numbers only. As an additional bonus, we go 
beyond the quantizations of the classical cases (3 = 1,2,4 and obtain continuous exponents 
(see Section IV for further discussion of this point). 

We provide real tridiagonal random matrix models for all /3-Gaussian (or Hermite) and /3- 
Wishart (or Laguerre) ensembles, and we discuss the possibility of constructing a real matrix 
model for the /3-MANOVA (or Jacobi) ensembles. 

We obtain our results by extrapolating the classical cases, thereby providing concrete 
models for what have previously been considered purely theoretical distributions. 

In Section II we establish results for symmetric tridiagonal matrices, and we use them to 
construct tridiagonal models for the 0- Hermite ensembles. Along the way, we obtain a short 
proof based on Random Matrix Theory for the Jacobian of the transformation T — * (q,X), 
where T is a symmetric tridiagonal matrix, A is its set of eigenvalues, and q is the first row 
of its eigenvector matrix. In Section III we construct tridiagonal models for the /3-Laguerre 
ensembles, by building on the same set of ideas that we use in Section II. In Section IV 
we present some immediate applications of the new classes of ensembles and we discuss the 
/?-Jacobi ensembles and other interesting open problems. 

We display our random matrix constructions in Table |l[ 



Hermite Matrix 
nGN 




/ N(0,2) X (n-l)P \ 
X(n-l)p N(0,2) X(n-2)/3 

X 2(3 N(0,2) xp 
\ X p N(0,2) J 


Laguerre Matrix 

m G N 
a e R 
a > §(m - 1) 


Lp = B„B 

Bp ~ 


J, where 

( X2a \ 
\ XP X2a~P{m-l) J 



Table 1: Random Matrix Constructions 
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2 The /3-Hermite (Gaussian) Ensembles 



2.1 Motivation: Tridiagonalizing the GOE, GUE, and GSE 

The joint distribution fp(\) of the eigenvalues for the GOE, GUE, and GSE is 

/ /3 (A) = 4|A(A)|^e-^» A ?, (7) 

where (3 = 1, 2, 4 [18]. Here the Vandermonde determinant notation A(A) stands for — 

Xj), and c? H is given by (|2|). 

We will prove in Section II. B that the tridiagonal /3-Hermite random matrix displayed in 
Table |l] has the joint eigenvalue p.d.f. given by general (5 in (|7[). For motivation, we will begin 
with a quick "back-door" proof for /3 = 1 by tridiagonalizing the GOE; then we will extend 
the result to the GUE and GSE. 

To illustrate the proof and help the reader follow it more easily, we have included the 
diagram of Figure ||. 

Theorem 2.1. If A is an n x n matrix from the GOE, then reduction of A to tridiagonal 
form shows that the matrix T from the 1-Hermite ensemble has joint eigenvalue p.d.f. given 
by (0j with (3 = 1. 

" J . Here a n is a standard Gaussian, a; is a vector of (n — 1) 

i.i.d. Gaussians of mean and variance 1/2, and B is an (n— 1) x (n— 1) matrix from the GOE; 
a n , x and B are all independent from each other. 

Let H be any (n — 1) x (n— 1) orthogonal matrix (depending only on x) such that 

Hx = [\\x\\ 2 ... 0] T = ||x|| 2 ei , 

where e\ = [1, 0, . . . , 0] T . Then clearly 

1 W a n x T \ ( 1 \ _ / a n \\x\\ 2 e\ 



H J \ x B J \0 H T J \ \\x\\ 2 ei HBH T 

Since A is from the GOE and H depends only on x, we can readily identify the distributions 
of a n , \\x\\2 and HBH T (these three quantities are clearly independent). The entry a n is 
unchanged and thus a standard normal with variance 1. Being the length of a multivariate 
Gaussian of mean and entry variance 1/2, ||x||2 has the distribution A=Xn—i- It is worth 

mentioning that the p.d.f. of \\x\\ 2 is given by r ,n-i , y n ~ 2 e~ y2 . Finally, by the orthogonal 

invariance of the GOE, HBH T is an (n—1) x (n— 1) matrix from the GOE. 
Proceeding by induction completes the tridiagonal construction. 

Because the only operations we perform on A are orthogonal similarity transformations, 
which do not affect the eigenvalues, the conclusion of the theorem follows. □ 

We recall that matrices from the GOE have the following properties: 

_ I yj a? 

Property 1. the joint eigenvalue density is c\ |A(A)|e » ' fl~8|| ; 

Property 2. the first row of the eigenvector matrix is distributed uniformly on the sphere, 
and it is independent of the eigenvalues. 
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The second property is an immediate consequence of the fact that the eigenvector matrix 
of a GOE matrix is independent from the eigenvalues Jl8| , (3.1.3) and (3.1.16), pages 55-58], 
and has the Haar (uniform) distribution because of the orthogonal invariance. 

The following corollary is easily established. 

Corollary 2.2. If T is a matrix from the 1-Hermite ensemble, with eigendecomposition T = 
QAQ T , then the first row q of the eigenvector matrix Q is independent of A, and is distributed 
uniformly on the sphere. 

Proof. If A = QiAQj and T = HAH T , then Q = HQ\. Since each one of the reflectors 
which form H has first row e±, multiplication by H does not affect the first row of Q\. The 
conclusion follows. □ 

Reduction to tridiagonal form is a familiar algorithm which solves the symmetric eigenvalue 
problem. The special "reflector" matrix H used in practice for a vector x = [x\,... , x n -i] T is 



H = 1 -2 



T 

UU 



U T U 



where u = x ± x\ e\. This special matrix H is known as the "Householder reflector" (see [11 
page 209]). 

The tridiagonal reduction algorithm can be applied to any real symmetric, complex her- 
mitian, or quaternion self-dual matrix; the resulting matrix is always a real, symmetric tridi- 
agonal. Using the algorithm similarly on a GUE or GSE matrix one gets the following 

Corollary 2.3. When (3 = 2,4, reduction to tridiagonal form of matrices from the GUE, re- 
spectively GSE, shows that the tridiagonal 2-Hermite, respectively 4-Hermite, random matrix 
has the distribution given by Note that (5 "counts" the number of independent Gaussians 
in each entry of the matrix. 

Remark 2.4. The observation that Numerical Linear Algebra algorithms may be performed 
statistically is not new; it may be found in the literature (see Trotter fi3J\J , Silverstein 
Edelman j^j). 

2.2 Tridiagonal Matrix Lemmas 

In this section we prove lemmas that will be used in our constructions in Sections II. C and 
III.B. 

Given a tridiagonal matrix T defined by the diagonal a = {a n , . . . , a%) and sub-diagonal 
b = (b n -i, ... , &i), with all bi positive, let T = QAQ T be the eigendecomposition of T as in 



Theorem 2.12. Let q be the first row of Q and A = diag (A). 



Lemma 2.5. Under the assumptions above, starting from q and A, one can uniquely recon- 
struct Q and T. 



Proof. This is a special case of the more general Theorem 7.2.1 in Parlett [24]. □ 
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Remark 2.6. It follows that, except for sets of measure 0, the map T — > (q, A) is a bijection 
from the set of tridiagonal matrices of size n with positive sub- diagonal, to the set of pairs 
(q,\), with q a unit norm n- dimensional vector of positive real entries, and A a strictly in- 
creasingly ordered sequence of n real numbers. Let the bisection's Jacobian be denoted by J 

('-{«$})■ 

Our next lemma establishes a formula for the Vandermonde determinant of the eigenvalues 
of a tridiagonal matrix. 

Lemma 2.7. The Vandermonde determinant for the ordered eigenvalues of a symmetric tridi- 
agonal matrix with positive sub-diagonal b = (6 n _i, ... , &i) is given by 

n-l 

n k 

A(A)=ij(*i-v ' 1 



t<3 n Q 

i=i 



n 



where (qi, . . . , q n ) is the first row of the eigenvector matrix. 

(k) 

Proof. Let A- , i = 1 ... k, be the eigenvalues of the k x k lower right-corner submatrix of T. 

Then Pk(x) = Yli=i( x ~ ^[ ) 1S the associated characteristic polynomial of that submatrix. 
For k = 1, . . . , n we have the three-term recurrence 

P k {x) = {x- a h )Pk-i(x) - bl_ x P k „ 2 {x) , (8) 

and the two-term relation 

n lAf^Af-^Hnip^cAf^nnip^Af- 1 ))! . & 
i < i < k i=i j=i 

1 < fk - 1 
From (||) we get 

i nW?- 1 ^ = n^- 2 (Af-^)i . (io) 

i=i i=i 

By repeatedly applying (||) and (|2.9|) we obtain 

n— 1 n—2 

u\p n (\t i] )\ = n i^-i(Ai n - 2) )i (ii) 
i=i i=i 

n-2 

(2(n-l) ,2(n-2) | TT p 



e 2 - 2) m^-3(Ar 2) )i (12) 
t=i 

(13) 

n-l 

n &r • (14) 



i=i 



s 



Finally, we use the following formula due to C.C. Paige, found in [24] as the more general 
Theorem 7.9.2: 



It follows that 



which proves the result. 



Pn-l(\) 













n 
n 

2 _ i=l 



(n)> 



n 

1=1 



A(A)S 



71-1 

n 



(15) 



(16) 
□ 



Remark 2.8. T/ie Vandermonde determinant formula of Lemma 2.1 can also be obtained 
from the Heine formula, as presented in Deift H/, page 44- 

The next lemma section computes the Jacobian J by relating the tridiagonal and diagonal 
forms of a GOE matrix, as in Figure |2[ 

A = QAQ T 

[GOE, Wishart real] 



Tridiagonal 
Reduction 




Eigenvalues & 
First Eigenvector Row 



Tridiagonal 
1-Hermite, 1-Laguerre 
ensemble 



(q,k) 



q, X as in 
Theorem 2.1, 3.1 



Figure 2: A dense symmetric matrix A can be tridiagonalized (left side) or diagonalized (right 
side). In brackets, we provide the distributions starting with that of A (GOE or Wishart real). 



Lemma 2.9. The Jacobian J can be written as 

n-l 

n bt 

j - i=i 

J — n 
i=l 
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Proof. To obtain the Jacobian, we will study the transformation from GOE to 1-Hermite 
ensemble (see Figure |2|). Note that J does not depend on (3; hence computing the Jacobian 
for this case is sufficient. 

Let T be a 1-Hermite matrix. We know from Section II. A that the eigenvalues of T are 
distributed as the eigenvalues of a symmetric GOE matrix A, from which T can be obtained 
via tridiagonal reduction (T = HAH T for some orthogonal H, which is the product of the 
consecutive reflections described in Section II. A). 

The joint element distribution for the matrix T is 

_I V a 2 n - T b 2 on-1 

H(a,b) = c a)b e i=1 [[^ e *=i , where c a>b 

i=l 



n-l 

(2vr)«/2 n r(i) 



Let 

da = A" =1 daj, d& = A^db i: dX = A" =1 Aj, 

and cig be the surface element of the n-dimensional sphere. Let fi(a(q, X),b(q, A)) be the 
expression for n(a,b) in the new variables q, X. We have that 

/J,(a,b) da db = J n(a(q,X),b(q,X)) dq dX = u(q,X)dqdX. (17) 

We combine Properties 1 and 2 of Section II. A to get the joint p.d.f. v(q, X) of the 
eigenvalues and first eigenvector row of a GOE matrix, and rewrite it as 

u(q, A) dq dX = n\ c\j j— ^— A(A) e i ' dq dX. 

We have introduced the n! and removed the absolute value from the Vandermonde, because 
the eigenvalues are ordered. We have also included the distribution of q (as mentioned in 
Property 2, it is uniform, but only on the all-positive 2 _n th of the sphere because of the 
condition > 0). 

n 

Since orthogonal transformations do not change the Frobenius norm ||j4||_f = ^2 a ij °f a 

1,3=1 



matrix A, from (17), it follows that 



J 



Kg,A) = n!ck-^ A(A) 
/j,(a, b) c a , b 



i=l 



i 



All constants cancel, and by Lemma 2.7 we obtain 



J 



n-l 

n bi 

i=l 
n 

n * 

8=1 



Note that we have not expressed /i(a, b) in terms of g and A in the above, and have thus 
obtained the expression for the Jacobian neither in the variables q and A, nor a and b, solely; 
but rather in a mixture of the two sets of variables. The reason for this is that of simplicity. □ 
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Remark 2.10. Our derivation of the Jacobian is a true Random Matrix derivation. Alter- 
nate derivations of the Jacobian can be obtained either via symplectic maps or through direct 
calculation. We thank Percy Deift and Peter Forrester, respectively, for having shown them 
to us. 

The last lemma of this section computes one more Jacobian, which will be needed in 
Section III.B. 

Let B be a bidiagonal matrix with positive diagonal x = (x m , ... ,x\) and positive sub- 
diagonal y = (y m -i, ■■■ ,yi)- Let T = BB T ; denote by a = (a m , ...ai) and b = (6 m _i, ... ,61) 
the diagonal, respectively the sub-diagonal part of T. Since T is a positive definite matrix, the 
transformation B — > T is a bijection from the set of bidiagonal matrices with positive entries 
to the set of positive definite tridiagonal matrices. 



Lemma 2.11. The Jacobian J, 



is 



m _j 

J(B^T) = (2 m ZiIJsi) • 

i=2 



Proof. We compute J^b^t) from the formula 

dx dy = J(b^t) da db , 

where dz = Aidzi for all z £ {a, b, x, y}. 
We have that 



ai 
h 



x r 



2 1 2 
Vi +Xi 

ViXi+i , 



(18) 
(19) 
(20) 



for alH = m — 1, m — 2, . . . ,1. 

Hence by computing differentials we get 

da m — 2x m dxin 

dai = 2(xj dxi + yi dyi), Mi = m — 1, m — 2, . . . ,1 
dbi = x i+ i dyi + yi dx i+ i, Vz = m - 1, m - 2, . . . , 

from which the formula follows. 



□ 



2.3 The Eigendistribution of the /3-Hermite ensemble 

Let #0 be a random real symmetric, tridiagonal matrix whose distribution we schematically 
depict as 



Ha 



1 

71 



( N(0,2) X (n-l)/3 

X(n-l)/3 #(0,2) X(n-2)/3 



\ 



X2/3 JV(0,2) X/3 

X/3 iV(0,2) / 
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By this we mean that the n diagonal elements and the n — 1 sub-diagonals are mutually 
independent, with standard normals on the diagonal, and -^Xk/3 on the sub-diagonal. 

Theorem 2.12. Let Hp = QAQ T be the eigendecomposition of Hp; fix the signs of the first 
row of Q to be non-negative and order the eigenvalues in increasing order on the diagonal 
of X = diag(A). Then X and q, the first row of Q, are independent. Furthermore, the joint 
density of the eigenvalues is 



uw=4 ni A *- A / e "" £Ai2 



4 \A(X)f e 



4E a? 



i<j 



and q = (qi, . . . , q n ) is distributed as (xp, ■ ■ ■ , Xp)i normalized to unit length. 



Proof of Theorem 2.H . 

Just as before, we denote by a = (a n , ... , ai) the diagonal of Hp, and by b = 
the sub-diagonal. The differentials da, db, dq, dX are the same as in Lemma [2.9| 
For general (3, we have that 



'M-i, 



n-l 



n-l 



(dHp) = n(a,b) da db = c afi f[ tff - 1 e"^ l|Tl||F da db = c U}b J JJ b\ fi ~ x e-^ T ^ F dq dX , 
where 



k=l 



k=l 



in— 1 



C a ,b 



With the help of Lemmas 2.7 and 2.9 this identity becomes 



(dHp) 



-r-rn-l 7 n-l 

C, h U$F^- II bf-'e-hW^ dq dX 



^a,b i-i-n I I 1 



-r-rn-l i k/3 n 

111=1% i=\ 



(21) 
(22) 



Thus 



(dHp) = ( rfflqt 1 dq) (n! <? H A(A) /3 e"^ ! A? dX 



i=l 



Since the joint density function of q and A separates, q and A are independent. Moreover, 
once we drop the ordering imposed on the eigenvalues, it follows that the joint eigenvalue 
density of Hp is (P H | A(A)|^e~ 2 ^ Xi , and q is distributed as (xp, ■ ■ • , X/?), normalized to unit 
length. From (p2[), it also follows that 



2 n - 1 r(fn) 

r(f) 



(23) 



□ 
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3 The /3-Laguerre (Wishart) Ensembles 



3.1 Motivation: Tridiagonalizing the Wishart ensembles 

The preceding section gives tridiagonal random matrix models for all /?-Hermite ensembles. In 
the following we define the /3-Laguerre ensembles, and give tridiagonal random matrix models 
for them. 

The Wishart ensembles have joint eigenvalue density 

m 

MX) = 4 a \A(X)f JjAT'e-^i V2 , (24) 

again with a = ^n, p = 1 + ^(m — 1), and with (3 = 1 for real, respectively (3 = 2 for complex. 
Here c^' a as the same as in (||). 

From now on p will always denote the quantity 1 + ^(m — 1), following the notation of 
Muirhead for (3 = 1 [21| (chapter 7) and Forrester (Forrester uses 1 + -(m — 1), where 



a = 2/(3 is the Jack parameter). Its presence is implicit in the p.d.f. of all /3-Laguerre 
ensembles; hence we will identify the ensembles by (3 and by a (we call the latter the "Laguerre" 
parameter, generalizing from the univariate case (3 = 1, m = 1). 

As in Section II. A, we will provide the most basic case for our construction: the case (3 = 1 
and Wishart real exponent n ~ r ^~ 1 (also referred as the case (3 = 1 and Laguerre parameter 
a = f). 

Theorem 3.1. Let G be an m x n matrix of i.i.d. standard Gaussians; then W = GG T 
is a Wishart real matrix. By reducing G to bidiagonal form B one obtains that the matrix 
T = BB T from the 1-Laguerre ensemble of Laguerre parameter a = § (defined as in Table ^) 
has the joint eigenvalue p.d.f. given by n2A). 



Proof. We write 

G 



T 

X 

Gx 



with x T a row multivariate standard Gaussian of length n and G\ a (m — 1) x n matrix of i.i.d. 
standard Gaussians. Let R be a right reflector corresponding to the vector x T (R T x = \ \x\ I2 ej) 
which is independent of G\. Hence G\R is a matrix of i.i.d. standard Gaussians. 

Write G\R = [y, G2], where y is a column multivariate standard Gaussian of length m — 1 
and G2 is a (m — 1) x (n — 1) matrix of i.i.d. standard Gaussians. Let L be a left reflector 
corresponding to y (Ly = \ \yW2 &i) which is independent of G2. Then we have that 

1 \ / ||x|| 2 

L J V 1Mb ei LG 2 

As we have seen before, ||x||2 is distributed like Xn-i, IMk is distributed like Xm-i, and LG2 
is a matrix of i.i.d. standard Gaussians (since L and G2 are independent). 
We proceed inductively to finish the bidiagonal construction of B. 

Because the operations we have performed on G are orthogonal left and right multiplica- 
tions, which do not affect the singular values, it follows that the singular values of G and B 
are the same. Since the squares of the singular values of G, respectively B, are the eigenvalues 
of W, respectively T, the conclusion of the theorem follows. □ 
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Remark 3.2. The bidiagonalization process presented above is part of a familiar Numerical 
Linear Algebra algorithm for computing the singular values of a matrix. 

Corollary 3.3. The same process of bidiagonalization performed on G,, a matrix of i.i.d. 
standard complex (standard quaternion) Gaussians, shows that the matrix W = GG T and the 
matrix T from the 2-Laguerre (A-Laguerre) ensemble of parameter a = n (a = 2n) has the 
joint eigenvalue p.d.f. given by (p^h In all three cases (real, complex, quaternion) we say that 
T represents the tridiagonalization of the Wishart (real, complex, quaternion) ensemble. 

In the next section we prove the general form of the theorem. 



3.2 The Eigendistribution of /3-Laguerre Ensemble 



Let 



Br, 



( X2a 

Xp(m-1) X2a-/3 

V 



\ 



Xp X2a-P(m-1) ) 



by this meaning that all of the 2m — 1 diagonal and subdiagonal elements are mutually 
independent with the corresponding \ distribution. 
Let Lp 



BpBj be the corresponding tridiagonal matrix. 



Theorem 3.4. Let Lp = QAQ T be the eigendecomposition of Lp; fix the signs of the first row 
of Q to be non-negative and order the eigenvalues increasingly on the diagonal of A. Then A 
and the first row q of Q are independent. Furthermore, the joint density of the eigenvalues is 



' =1 Ai/2 



1=1 



where p = 1 + ^ (m — 1), and q is distributed as (xp-i ■ ■ ■ > Xp) normalized to unit length. 



Proof of Theorem [O. We will use throughout the results of Lemma 2.7, Lemma 2.9 



Lemma 2.1l| and Remark 2J3, which are true in the context of tridiagonal symmetric matrices 
with positive sub-diagonal entries. By definition, Lp is such a matrix. 

We will again use the notations of Lemma 2J} and 2.11 for the differentials da, db, dq, dX, 
dx, and dy. 

We define (dBp to be the joint element distribution on Bp 

m—l rn—l 

(dBp) ee vi{x,y)dxdy = W^t 1 ^' 2 Wv^e^' 2 dx dy . 

i=0 i=l 

By using Lemma [24| we obtain the joint element distribution on Lp as 

(dLp) ee J B ^ T fi(x,y) dx dy (25) 

m— 2 m—l 
= 2-™C Xiy x \^m-l)-2 e - x y 2 jj ^-3^/2 J] y Pi -l e -,|/2 ^ dy (26) 



i=0 



i=l 
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where 



We rewrite (26) in terms of x, y, A, and q: 

1 li=l Qi 

m— 2 m— 1 

x n ^ (m - i) - 3 n ^r 1 * ^ 

i=0 t=l 

I U=l 1i 

Since the Vandermonde with respect to b and q and the ordered eigenvalues A can be 
written as 



A(A) - lli=1 1 



lli=l?< 

it follows that 

lir:.'( 

A(A) 



n*=i % 

This means that we can rewrite 

Pi 



lli=l?8 i=i j=0 



777 — 1 



™-i ™ 1 /mi \ 2a— fl(m— 1)— 2 

-2 /n v 1 ..2/" m-1 /m-1 \ <"v ; 



2- m c i2/ e" & £ * /2 A(A)^ J] rf -1 ( II ) * dA 



8=1 \ 2=0 



The trace and the determinant are invariant under orthogonal similarity transformations, 
so tr(Lg) = tr(A), and det(Lg) = det(A). This is equivalent to 



m— 1 



8=0 



m— 1 


m 

= 


8=1 


8 = 1 


m— 1 

TT x 2 

J J_ -^m— j 

8=0 


m 

= n*< 

8=1 



Using this, and substituting p for 1 + f (m — 1), we obtain that 
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where Cq is the same as in (|2 



From the above we see that q and A are independent, and once we drop the ordering 
the joint eigenvalue density is given by the /3-Laguerre ensemble of parameter a, while q is 
distributed like a normalized vector of X/3' s - 



This concludes the proof of Theorem 3.4 



□ 



4 Applications and Open Problems 

As we mentioned in Section |l], we believe that there should be many applications for the new 
tridiagonal ensembles. We illustrate here some (in Section IV. A), in the hope that researchers 
will find many more. Some of the applications we believe are new results (Applications 1, 3, 
5, and 6), and some are simplifications of known results (Applications 2 and 4). 

We discuss the open problem of constructing a matrix model for the /3-Jacobi ensembles 
in the beginning of Section IV. B. To facilitate the finding of new results, we conclude with a 
few open "general /3-ensemble" problems. 



4.1 Applications 

1. Interpolating Laguerre exponents. Our /3-Laguerre ensembles have "continuous" Laguerre 
parameters a which, even in the cases (3 = 1,2,4, interpolate the Wishart parameters. 

Though /3-Laguerre ensembles with general ( "continuous" ) parameter a have been stud- 
ied by many researchers (Q, flOU , |2l|| ), no non-quantized matrix realizations (i.e. explicit 
random matrix models) of /3-Laguerre ensembles are found in the literature. 
By "quantized" we mean that the exponent a is either an even integer, an integer, or a 
half- integer (depending on the value of (3). In particular, all models corresponding to a 



Laguerre (or Jacobi) weight found in [g3| and |T^] are quantized. 

Thus, our /3-Laguerre random matrix constructions extend the pre-existing ones in two 
ways: through (3 and through the Laguerre parameter a. 

The expected characteristic polynomial. The result below might be seen as an extension 



of the classical Heine theorem (see Szego [25] and Deift |q]) which has (3 = 2. Note that 



for (3 7^ 2, A(A)^ can no longer be written as the determinant of a Vandermonde matrix 
times its transpose, and the proof cannot be duplicated. 

The same result is found in a slightly more general form in ||, and its Jacobi case was 
first derived by Aomoto [jij. 

Theorem 4.1. The expected characteristic polynomial P n (y) = det(yl n — S) over S in 
the (3-Hermite, respectively (3-Laguerre, ensembles are proportional to 

( V \ it~ n ( V 

Hn J ' respectively to L£ 

f -n 

Here H n and Ln are the Hermite, respectively Laguerre, polynomials, and the constant 
of proportionality accounts for the fact that P n (y) is monic. 
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Proof. Both formulas follow immediately from the 3-term recurrence for the character- 
istic polynomial of a tridiagonal matrix (see formula (||)) and from the independence of 
the variables involved in the recurrence. □ 

3. Expected values of symmetric polynomials. Using the three-term recurrence for the char- 



acteristic polynomial of a tridiagonal matrix, we obtain Theorem 4.2. 

Theorem 4.2. Letp be any fixed (independent of (3) multivariate symmetric polynomial 
on n variables. Then the expected value of p over the (3-Hermite or (3-Laguerre ensembles 
is a polynomial in (3. 

We remark that it is difficult to see this from the eigenvalue density. 
Proof. The elementary symmetric functions 

6j (xi , X2 j • • • j X n ") — ^ %ji Xj 2 • • • Xj i i — 0,1,... , Tl , 

i<ii<... <n<n 

can be used to generate any symmetric polynomial of degree n (in particular p) . 

The ei evaluated at the eigenvalues of a matrix are the coefficients of its characteristic 
polynomial, and hence they can be written in terms of the matrix entries. Thus p can 
be written as a polynomial of the nx n tridiagonal matrix entries (which corresponds to 
the Hermite, respectively Laguerre cases). 

To obtain the expected value of p over the /3-Hermite or /3-Laguerre ensemble, one can 
write p in terms of the corresponding matrix entries, use the symmetry to condense the 
expression, then replace the powers of the matrix entries by their expected values. 

The diagonal matrix entries are either normal random variables in the Hermite case or 
sums of x 2 random variables in the Laguerre case. The subdiagonal entries appear only 
raised at even powers in the et and hence in p (this is an immediate consequence of the 
three-term recurrence for the characteristic polynomial, (||)). Since all even moments of 
the involved x distributions are polynomials in /3/2, it follows that the expectation of p 
will be a polynomial in /3. □ 

As an easy consequence we have the following corollary. 

Corollary 4.3. All moments of the determinant of a (3-Hermite matrix are integer- 
coefficient polynomials in (3/2. 

Proof. Note that even moments of the %Bi distribution are integer-coefficient polynomials 
in (3/2, and that the determinant is e n . □ 

4. A new proof for Hermite and Laguerre forms of the Selberg Integral. Here is a quick proof 
for the Hermite and Laguerre forms of the Selberg Integral (|l8|]), using the /3-Hermite, 
respectively, /3-Laguerre ensembles. 
The Hermite Selberg integral is 



I H ((3,n)= f \A(\)fe-^ x * /2 d\ 
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We have that 



I H ((3,n) = n\ ( I A(A) /3 e - E -i A ? /2 dX ) 

\Jo<Ai<...<A n <oo / 




where Cq is as in (23). We introduce the n! because in the first integral we have ordered 
the eigenvalues; signifies that all qi are positive. 

Note that Cg can easily be computed independently of the /3-Hermite ensembles. 

Using the formula for the Vandermonde given by Lemma 2/7 , the formula for the Jacobian 
J given in Lemma 2.9, and the fact that the Frobenius norm of a matrix in the tridiagonal 
1-Hermite ensemble is the same as the Frobenius norm of its eigenvalue matrix, one 
obtains 

1*08, n) = n\£\ [U =\ Qi LUi K TT^-ie-Etx/g-SEW/* da db 

q Vx ( o,oo)^ n^i 1 h nr =1 «f tr 

n-1 



n\ c^(2^)"/ 2 TT / b^e-'i db 

7=1 J {0,°o) 



n 



(r(D) fi 2 4 



The same reasoning yields the Laguerre Selberg Integral formula 



I 



f3,a,n 
L 



1 

C L 



5. Moments of the discriminant. The discriminant of a polynomial equation of order m 
is the square of the Vandermonde determinant of the m zeroes of the equation. Thus, 
the discriminant of the characteristic polynomial of a /?-Hermite of /3-Laguerre ensemble 
matrix is simply D(\) = A (A) 2 . 

A simple calculation shows that the kth. moment of D(\) is 



]J "7 Vr^" ' respectively, 

2 fcm(m-l) V -rfj'Jfci (o - f(m -./ i I/.,,- i , 



" l (l + fi)fci ( a-f(m -j)) 

(1 + f )* 



/3+2fc, o+fe(m-l) 11 /i , g 



where n and m are the matrix sizes for the Hermite, respectively, Laguerre cases, and 
the rising factorial (x)k = T(x + k)/T(x). 

Using the Selberg integral, one obtains that the moments of the discriminant for the 
/3-Jacobi case are 

c? ai '° 2 _ (1 + f j)fej (ai-f(m-j)) fc(3 --i) M(m-j)fa 

^ +2fc , a 1+ Mm-l), a2 +fc (m- 1 ) 11 (1 + j )jfc (oi + a2 _ | (m _ , ))fe(m+j ._ 2) 
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6. Software for Application 2. Computing eigenvalue statistics for the (5-ensembles. Appli- 
cation 2 suggests that integrals of the form 



may be evaluated with software. 

One example of this would be computing moments of the determinant over the /3-Hermite 



ensemble. There are explicit formulas for the cases (3 = 1,2 and 4, due to Mehta [19] 
and to Delannay and Le Caer ||, which can be used to evaluate these moments. 

In the absence of a closed-form, explicit formula, like the one for (3 = 1 provided in 
||, the computation of these moments cannot be made polynomial; thus it is inherently 
slow. 

For the general f3 case, one can compute the moments in terms of a multivariate Hermite 
polynomial evaluated at (see Q, |2|). Using this technique, the complexity of the 
computation exceeds by far that of symbolically taking the determinant of a tridiagonal 
matrix, expanding the power, and replacing all powers of the entries by their expected 
values (which are all known). Writing a Mathematica code to implement this algorithm 
is an easy exercise, and such a code would allow the author to compute these moments 
in a reasonable amount of time, provided that the product between the power and the 
size of the matrix is not very large. A template for a special case when (3 = 1 can be 
found in II Appendix A] . 



4.2 Open Problems 

1. [3-Jacobi (MANOVA) Ensembles. Sections II and III of the paper provide tridiagonal 
matrix models for the /3-Hermite and /3-Laguerre ensembles. The natural question is 
whether such models exist for the last member of the classical triplet, Jacobi. The (3- 
Jacobi ensembles have been intensively studied as theoretical distributions, especially in 
connection with Selberg-type integrals and Jack (or Jack-Selberg) polynomials (see Q, 
fl5| , [pi , H). Finding a random matrix model that corresponds to them would be of 
much interest. 

If the two matrix factorizations problems that are associated with the Hermite and 
Laguerre ensembles are the EIG and the SVD, the one associated with the Jacobi should 
be the QZ (the generalized symmetric eigenvalue problem). This idea is supported by 
the fact that the MANOVA real and complex distributions, which correspond to the 
Jacobi /3 = 1,2 ensembles, are indeed connected to the QZ algorithm. A good reference 
for QZ is 0. 

Though we have not studied this problem sufficiently, we believe that a concrete (per- 
haps sparse, perhaps tridiagonal) matrix model may be constructed for the /?-Jacobi 
ensembles. 

2. Level densities. The level density of an ensemble is the distribution of a random eigen- 
value of that ensemble (and by the Wigner semicircular law we know that the limiting 
distribution as n — > oo of such an eigenvalue is semicircular). The three functions found 
to be the level densities of the Gaussian models depend on the univariate Hermite poly- 
nomials. 
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Recently, Forrester [10] has found a formula for the level densities of the /3-Hermite 
ensembles which works for (5 an even integer. This formula depends on a multivariate 
Hermite polynomial. 

Finding a unified formula for the general (3 case would be of interest. 

Level spacings. The level spacings are the distances between the eigenvalues of an en- 
semble, usually normalized so that the average consecutive spacing is 1. These spacings 
have been well-studied in the case of the Gaussian ensembles {(3 = 1, 2,4). The limiting 
probability density of a random spacing in these cases is known in terms of spheroidal 



functions (see [18]). 

A surprising connection exists between the limiting probability density of a GUE random 
spacing and the probability density of the zeroes of the Riemann zeta function. Inspired 
by the theoretical work of Montgomery ( |2(| ) , Odlyzko ( |22|] ) has shown experimentally 
that the two probability densities are very close; the subsequent conjecture that the two 
probability densities coincide has been named the Montgomery-Odlyzko law. 

To the best of our knowledge, the level spacing of the general /3-Hermite ensembles has 
not been investigated. 

4. Bulk and edge scaling limits. Finally, a very important application would be the gener- 
alization of the bulk and edge scaling limits for the GOE, GUE and GSE obtained by 
Tracy and Widom (the latter are known as the Tracy- Widom distributions F\, F2 and 
F 4 ). 

The edge scaling limit refers to the distribution of the largest eigenvalue of a matrix in the 
ensemble; the bulk scaling limit refers to the distribution of an eigenvalue in the "bulk" 



of the spectrum. For a reference, see [29], [30] or The Tracy- Widom distributions 



are defined in terms of Painleve functions, which are solutions to certain differential 
equations, with asymptotics given by Airy functions. For a good treatment of Painleve 
equations in relationship with Gaussian (Hermite), Laguerre, and Jacobi random matrix 
models, see Pierre van Moerbeke's notes (|32"| , Section 4]. Recently, Johnstone [|l4|] has 
found that the limiting distributions F\ and F2 apply to real (respectively complex) 
Wishart matrices. 
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